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We present a relativistic model describing a thin disk surrounded by a halo in presence of an 
electromagnetic field. The model is obtained by solving the Einstein-Maxwell equations on a partic- 
ular conformastatic spacetime background and by using the distributional approach for the energy- 
momentum tensor. A class of solutions is obtained in which the gravitational and electromagnetic 
" potentials are completely determined by a harmonic function only. A particular solution is given 

that is asymptotically flat and singularity-free, and satisfies all the energy conditions. 
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I. INTRODUCTION 

A large number of galaxies and other astrophysical systems have extended mass distributions surrounded by a 
material halo. For practical reasons one can assume that many of these systems preserve axial symmetry and, 
therefore, they can be modeled in terms of relativistic thin disks with exterior halos, in particular when the gravity 
field is strong enough. Disks may also be used to model accretion disks, galaxies in thermodynamic equilibrium and 
the superposition of a black hole and a galaxy. Disk sources for stationary axially symmetric spacetimes with magnetic 
fields arc also of astrophysical importance mainly in the study of neutron stars, white dwarfs and galaxy formation. 
To describe the gravitational and electromagnetic fields of such configurations, we will use general relativity and 
Maxwell's theory. Consequently, we are interested in deriving and analyzing exact solutions of the Einstein-Maxwell 
equations. On the other hand, the study of axially symmetric solutions of the Einstein and Einstein-Maxwell field 
equations corresponding to shells and disk-like configurations of matter, apart from its astrophysical relevance, has a 
clear purely mathematical interest. 

Exact solutions that have relativistic static thin disks as their sources were first studied by Bonnor and Sackfield 
[l[ and Morgan and Morgan 0, Q ■ Subsequently, several classes of exact solutions corresponding to static [i|-fi3| and 
stationary [l44f 81 ] thin disks have been obtained b y d ifferent authors. The superposition of a static or stationary thin 
disk with a black hole has been considered in [L9l - [27j . Thin disks around static black holes in a magnetic field have 
been studied in 1281. Relativistic disks embedded in an expanding Friedman-Lemaitre-Roberton- Walker universe have 
been studied in [29|. and perfect fluid disks with halos in [30j . Furthermore, the stability of thin disks models has 
been investigated using a first order perturbation of the energy- momentum tensor in [3l| . On the other hand, thin 
disks have been discussed as sources for Kerr-Newman fields [12, , magnetostatic axisymmetric fields [U, HH , and 
conformastatic and conformastationary metrics j36l - l38| . Also, models of electrovacuum static counterrotating dust 
disks were presented in [39[ , charged perfect fluid disks were studied in |4(j , and charged perfect fluid disks as sources 
of static and Taub-NUT-type spacetimes in [4f|, [42j. Also, monopole and dipole layers in curved spacetimes were 
analyzed in p3| , and electromagnetic sources distributed on shells in a Schwarzschild background in [44| . 

Now, the thin disks with magnetic fields presented in [32l - l34| were obtained by means of the well-known "displace, 
cut and reflect" method that introduces a discontinuity in the first-order derivative of an otherwise smooth solution. 
The result is a solution with a singularity of the delta-function type in the entire z = hypcrsurface, and so it can be 
interpreted as an infinite thin disk. On the other hand, solutions that can be interpreted as thin disks of finite extension 
can be obtained if an appropriate coordinate system is introduced. A coordinate system that adapts naturally to a 
finite source and presents the required discontinuous behavior is given by the oblate spheroidal coordinates. Some 
examples of finite thin disks obtained from vacuum solutions expressed in these coordinates can be found in references 
[J S 0; 1 an d from electrovacuum solutions in reference [38| . 

In a previous work [38| , we presented an infinite family of conformastatic axially symmetric charged dust disks of 
finite extension with well-behaved surface energy and charge densities. These disks have a charge density that is equal, 
up to a sign, to their energy density, and so they are examples of the commonly named "electrically counterpoised 
dust" equilibrium configuration. The energy density of the disks is everywhere positive and well-behaved, vanishing 
at the edge. Furthermore, since the energy density of the disks is everywhere positive and the disks are made of dust, 
all the models are in a complete agreement with all the energy conditions, a fact of particular relevance in the study 
of relativistic thin disks models. In the present paper, we extend these studies to obtain a model corresponding to 
a system that is composed of a thin disk and an exterior halo. The main purpose of this work is, then, to extend 
the previous electric field to include an electromagnetic field, and the previous "isolated" thin disk to include a thin 
disk-halo system. 

In this work, we present a relativistic model describing a thin disk surrounded by a halo in presence of an electro- 
magnetic field in a conformastatic spacetime (45l - l47j . Note that we take the definition in [46j as standard, following the 
original terminology by Synge (45j : conformastationary are those stationary spacetimes with a conformally flat space 
of orbits and the conformastatic comprise the static subset. The model is obtained by solving the Einstein-Maxwell 
equations and by using the distributional approach under the assumption that the energy-momentum tensor can be 
expressed as the sum of two distributional contributions, one due to the electromagnetic part and the other one due 
to a "material" part. In this way, explicit expressions for the energy, pressure, electric current and electromagnetic 
field are obtained for the disk region and for the halo region. In order to obtain the solutions, an auxiliary function 
is introduced that determines the functional dependence of the metric and the electromagnetic potential. It is also 
assumed that the auxiliary function depends explicitly on an additional function which is taken as a solution of the 
Laplace equation. A simple thin disk-halo model is obtained from the Kuzmin solutions of the Laplace equation. The 
energy-momentum tensor of the system agrees with all the energy conditions. 

The plan of our paper is as follows. First, in Section HU the conformastatic line element is considered. The 
procedure to obtain electromagnetostatic, axially symmetric, relativistic thin disks surrounded by a material halo is 
also summarized in this section. Section IIIII introduces a functional relationship dependence between the metric and 
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electromagnetic potentials and an auxiliary function in order to obtain a family of solutions of the Einstein-Maxwell 
equations in terms of convenient solutions of the Laplace equations, modeling relativistic thin disk-halo systems. Next, 
the eigenvalue problem for disk with halos is studied and a particular model of a disk with halo is obtained from the 
Kuzmin solutions of the Laplace equation. In Scction flll C[ the quadratic Riemann invariants and the electromagnetic 
invariants are studied, and the behavior of the Kuzmin disk with halo is analyzed. Finally, Section llVl is devoted to 
a discussion of the results. 



II. THE EINSTEIN-MAXWELL EQUATIONS AND THE THIN-DISK-HALO SYSTEM 

In order to formulate the Einstein-Maxwell equations for conformastatic axially symmetric spacetimes corresponding 
to an electromagnetized system, constituted by a thin disk and a halo surrounding the exterior of the disk, we first 
introduce coordinates x a = (t, ip, r, z) in which the metric tensor and the electromagnetic potential only depend on 
r and z. We assume that these coordinates are quasicylindrical in the sense that the coordinate r vanishes on the 
axis of symmetry and, for fixed z, increases monotonically to infinity, while the coordinate z, for r fixed, increases 
monotonically in the interval (—00,00). The azimuthal angle ip ranges in the interval [0, 2n), as usual 0, 0- We 
assume that there exists an infinitesimally thin disk, located at the hypersurface z = 0, so that the components of the 
metric tensor g ab and the components of the electromagnetic potential A a arc symmetrical functions of z and their 
first derivatives have a finite discontinuity at z — 0, accordingly, 

g a b(r,z)=g ab (r,-z), A a (r, z) = A a (r, -z) (1) 

in such a way that, for z 7^ 0, 

9ab,z( r > z ) = -9abA r > A a,z(r, z) = -A a ^ z (r, -z). (2) 

The metric tensor and the electromagnetic potential are continuous at z = 0, 

[gab] = 9ab\ z=0+ - 9ab\ z=Q - = 0, (3a) 

[A a ] = A a l =o+ -A a l =o _=0, (3b) 

whereas the discontinuity in the derivatives of the metric tensor and the electromagnetic potential can be written, 
respectively, as 

Jab=[9ab,z], (4a) 
Ca = [A a , z ], (4b) 

where the reflection symmetry with respect to z — has been used. Then, by using the distributional approach 
psT - llol ] or the junction conditions on the extrinsic curvature of thin shell [5ll - l53l ] , we can write the metric and the 
electromagnetic potential as 

9ab=9i b 0(z)+g- b {l-9(z)}, (5a) 
A a = A+9(z) + A-{l-9(z)}, (5b) 

and thus the Ricci tensor reads 

R ab = R+ b 0(z) + R- b {l - 6(z)} + H ab 5(z), (6) 

where 6{z) and 8{z) are, respectively, the Heaveside and Dirac distributions with support on z = 0. Here g^ b and R ab 
are the metric tensors and the Ricci tensors of the z > and z < regions, respectively, and 

Hab = \{l z Jl + l z b K - - g* z lab }, (7) 

where all the quantities are evaluated at z — + . In agreement with Eq. (|6]) , the energy- momentum tensor and the 
electric current density can be expressed as 

Tab = T+6{z) + T~ b {l - 6(z)} + Q ab S(z), (8a) 
J a = J+e(z) + J- {l - 6(z)} + I a 8(z), (8b) 
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where and J^ b are the energy-momentum tensors and electric current density of the z > and z < regions, 
respectively. Moreover, Q ab and I a represent the part of the energy-momentum tensor and the electric current density 
corresponding to the disk-like source. 

To describe the physical properties of an electromagnetized system constituted by a thin disk surrounded by an 
exterior halo, T^ b in (|8a[) can be written as 

T%, = B& + M±, (9) 
where E^ b is the electomagnetic energy-momentum tensor 

E ab = F ac F b c -^g ab F cd F cd , (10) 

with F ab = A b a — A a b and Mri is an unknown "material" energy-momentum tensor (MEMT) to be obtained. 
Accordingly, the Einstein-Maxwell equations, in gcometrized units such that c = 8itG = /-Iq = eq = I, arc equivalent 
to the system of equations 



Gab = Rab 0-9^ = Kb + M±, (11) 



Hab - ^9abH — Qab, (12) 

Ff, b = J£, (13) 
[F ah ]n b = r , (14) 

where H = g ah H a b and a = ^J—ga, being g the determinant of the metric tensor. Here, "[ ]" in the expressions [F ab ] 
denotes the jump of F ab across of the surface z = and n b denotes an unitary vector in the direction normal to it. 

Now, in order to obtain explicit forms for the Einstein-Maxwell equations corresponding to the electromagnetized 
disk-halo system, we take the metric tensor as given by the conformastatic line element (45l . |46| 

d,s 2 = -e 2<#, di 2 + e 2tp [r 2 dip 2 + dr 2 + dz 2 ], (15) 

where the metric functions <f> and ip depend only on r and z, and the electromagnetic potential is 

A a = (Ao, A, 0, 0). (16) 

We also assume that the electric potential Aq and the magnetic potential A are independent of t. For later use, 
the corresponding field equations and the components of the energy-momentum tensor are given explicitly in the 
Appendix. 

Furthermore, from Eq. (|13p we have for the current density in these regions, 

jo = _ re ^ {V 2 ^ - V(0 - VA }, (17a) 

Ji = -r" V-^-fV 2 ^ + V(0 - V)- VA - -A r }, (17b) 

r ' 

where, as we know, = J±y/—g. The "true" surface energy- momentum tensor of the disk, S a b, can be obtained 
through the relation 



Sab = J Qab 5{z) ds n = ^fg~T z Q a b, (18) 

where ds n = ^g Z z dz is the "physical measure" of length in the direction normal to the z — plane. Accordingly, for 
the metric t|l 5|) , the nonzero components of S a b are given by 

= 4e-^V, z , (19a) 

S 1 1 = 2e-^(0 + V) ) „ (19b) 

S\ = 2e-^ + V),z, (19c) 

where all the quantities are evaluated at z = + . The "true" current density on the surface of the disk, J a , can be 
obtained through the relation 

J a = j I a 5{z)^g— Z dz, (20) 
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or explicitly 

J°=e-^ +2 ^[A 0>z ], (21a) 
J 1 = -r- a e-»+[A t ,], (21b) 

where all the quantities are evaluated on the surface of the disk, and ^4 a ,z] denotes the jump of the derivative of A a 
across the surface 2 = 0. Now, in order to analyze the physical characteristics of the system it is convenient to express 
the energy-momentum tensor T ao and the electric current density in terms of an orthonormal tetrad. We will use the 
tetrad of the "locally static observers" (LSO) H3], i.e., observers at rest with respect to infinity, which is given by 



e 



(b) 



where 



{V a 7 W a ,X a ,Y a }, (22) 



(23a) 
(23b) 
(23c) 
(23d) 



ya 


= e (0) 


= e~ 




w a 


= e d) 


= r~ 






- p a 
~ e (2) 


= e~ 




ya 


- p a 
~ e (3) 


= e~ 





In terms of this tetrad M^ b and can be expressed as 



Mf = M± m V a V b + Mf l){1) W a W b + Mf 2)[2) X«X b + M$ ){3) Y a Y b (24a 
- M± )(1) {V a W b + W a V b ) + M± )(3) {X a Y b + Y a X b ) , 

T± ya , j± 



J a ± = -J±V a + J±W a . (24b) 



The explicit expressions for these quantities are given in the Appendix. 

In the same way, by using the LSO tetrad, the surface energy-momentum tensor and the surface current density of 
the disk as well as the electric current on the disk can be written in the canonical form as 

S ab = S mo) V a V b + S (l)(1) W a W b + S {2)(2) X a X b , (25a) 
J a = -J {0 )V a + J (1) W a , (25b) 



with 



5(o)(o) = -4e"^, z , (26a) 
5 (1)(1) = 2e"*(0 + V),, = 5(2X2), (26b) 
J(0) =-e-^[A , z ], (26c) 



J {1) =-r- l e-^[A,z], (26d) 



where we have used (|19[) and (|2T))) and again all the quantities are evaluated at z = 



+ 



III. THIN DISK WITH AN ELECTROMAGNETIZED MATERIAL HALO 



In the precedent section, we discussed a generalized formalism in which conformastatic axially symmetric solutions 
of the Einstein-Maxwell can be interpreted in terms of a thin disk placed at the surface z = surrounded by a 
distribution of electrically charged matter located in the z > and z < regions, whose physical properties can be 
studied by analyzing the behavior of 5(a) (fc), 3(a), ^u)(b) anc ^ ^uy T° this end, it is necessary to "choose" a convenient 
explicit form for the metric. It turns out that the assumption = —<f> leads to a considerable simplification of the 
problem. Indeed, the equations (|A6|) reduce to 
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M (o)(o) = Jr'WV 2 / - 5V/ • V/ - 2/VAq • VA - 2 r - 2 fVA ■ VA}, (27a) 
Mf Q){l) =-r- 1 /VA -Vi, (27b) 

M W(D = I- rl { y/ ' V/ ~ 2 ^ A ° ■ VA ° - 2r ~ 2 fVA ■ VA}, (27c) 
M% (2) = \r l {-{fl - f%) + 2/«, - Al z ) - 2r-if{A\ - A%)}, (27d) 

M 

M (3)(3) = - M tmv ( 27c ) 

Mg )(3) = 1/ : /,/ : + 4>,rA>,« - r- 2 f 2 A, r A iZ , (27f) 

■^(0) = ^/ 1/2v ' (/~ 'VA)), (27g) 
Jj) =-r 2 /- 1 / 2 V-(r- 2 /VA). (27h) 

In the same way, for the nonzero components of the energy-momentum and the current density on the surface of 
the disk (|26|) we have, respectively, 

S (0) (o) = 4(/ 1/2 ), z , (28) 

and 

J { o) = -[Ao,m], (29a) 

J(i) = -r-VlAJ. ( 29b ) 

where "[ ]" denotes the jump across of the disk, / = e 2 ^ and all the quantities are evaluated on the surface of the 
disk. We will suppose that there is no electric current in the halo, i. e., we assume that = 0. Then, if e v is a unit 
vector in the azimuthal direction and A is any reasonable function independent of the azimuth, one has the identity 

V • (r-% x VA) = 0. (30) 

Equation (|27h[) may be regarded as the integrability condition for the existence of the function A defined by 

r~ 2 /VA = r~% x VA, (31) 

or, equivalently 

-/- 1 VA = r"%xVA (32) 

Hence, the identity (f3T)| implies the equation 

V • (r X VA) = (33) 

for the new "potential" A(r, z). In order to obtain an explicit form of the metric and electromagnetic potential, we 
suppose that /, A and A depend explicitly on A. Then, we obtain from Q27gD 

(-r 7'A, + OVA • VA + A' Q W 2 X = 0, (34) 

where ()' denotes derivatives respect to A. Whereas, from (f33|) . we have 

-/~7'VA- VA + V 2 A = 0. (35) 

Consequently, by inserting ([33]) into (fM]) . we obtain Aq = 0, whose general solution is A = kiX + k 2 , where fei and fc 2 
are constants. We now proceed to determine the function A(r, z). Let us assume in (|35[) the very useful simplification 
f'f 1 = k, where k is an arbitrary constant. Then, / = k^e kx , and 

V 2 A = fcVA-VA, (36) 

where fc 3 is a constant. Furthermore, if we now assume the existence of a function U = k^e^^ + k 5 , with fc 4 and k 5 
being arbitrary constants, then 

V 2 [/ = -fcfc 4 e" feA (V 2 A-fcVA-VA) =0 (37) 
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and, consequently, A can be represented in terms of solutions of the Laplace equation: 



U - k 5 



V 2 U = 0. 



(38) 



Then, in order to have an asymptotically flat spacetime at infinity, we will only consider functions U that vanish at 
infinity and we must take fes = —k^k^ in Eq. (|38l) . 

On the other hand, from Eg. pip wc obtain the following relationship between A and A: 

\7A = A. r e r + A, z e z = r/" 1 ^ x (A, r e r + \ z e z ), (39) 

that is, A jr = — r/ _1 A iZ , and A tZ = rf \ r , or, in terms of U, A iT = k^rU^, and A_ z = —k^rU^, where kg = 
l/(kkzki). Then, using this solution for the nonzero components of M^ a ^ b y we have from (|27)) : 



± ul + u 2 z t 

M (o)(o) = 4([/ _ fc5 ; 2 { 3 / m, 


(40a) 


U 2 + U 2 

M (V)= 4(c ;-4 2{/ - M ' 


(40b) 


, u 2 r - u% 


(40c) 




(40d) 


— — M ± 

m (3)(3) - 1U (2)(2)- 


(40e) 



Furthermore, for the nonzero components of the energy-momentum tensor and the current density on the surface of 
the disk we have from (|28|) and ([29]) respectively, 



S( )(o) = 2(fc 3 fc 4 ) 1/2 



U, 



(u - k 5 yn 



and 



J(o) 
J(i) 



fe 



U-ks 



k(u-k 5 y 



with 



/ = 



fc 3 fc 4 



and 



2(kf + 1) 
fc 2 ' 



(41) 

(42a) 
(42b) 

(43) 



As we can see, 5(o)(o) is the only nonzero component of the surface energy-momentum tensor; we can then interpret it 
as the surface energy density of the disk e(r) as seen from any LSO. So, in order that the surface energy-momentum 
tensor will be in agreement with the energy conditions, we only need to require that U lZ \n+ > and take k^k^ > 0. 
Likewise, we interpret a(r) = Jr \ and I = Jru as the charge density and electric current density on the surface of 
the disk, respectively. 



A. The eigenvalue problem for the energy-momentum tensor of the halo 

When the tensor M^^) m the LSO tetrad is diagonal, its interpretation is immediate. In our case, however, 
^(2) (3) ^ ^ an d, therefore, it is necessary to rewrite M^\,^ in the canonical form. To this end, we must solve the 
eigenvalue problem for M^^, i.e., 

Mf k){b) e\=\ A%am e ) A , (44) 
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and express the physically relevant quantities in terms of the eigenvalues and eigenvectors. The solution of the 
eigenvalue problem in the LSO orthonormal tetrad leads to the eigenvalues 



A " - ~ A/ (o)(o)' 
Ai = M± (1) 

A± = ±VD, 



and the corresponding eigenvectors are given by 



tO) 
s 


= v (a) 


= (1,0,0,0), 


&\ 


= lW 


= (0,1,0,0), 


t (a) + 


_ y(a) 


= #(0,0,1,-0;), 




= z (a) 


= JV(0, 0,^,1), 



(45a) 
(45b) 
(45c) 



(46a) 
(46b) 
(46c) 
(46d) 



where 



D = 



- ( M (2)(2)) + i M (2) (3)) ' 



M, 



to 



N 



(2)(2) 



M, 



(2)(3) 



In terms of the tetrad £ (q) a = {^ a U W ^ ( "U W }, the tensor M,± )(f)) can be written in the canonical form 



M 



(a)(6) 



£^(a)V( 6 ) +PlX( a )X (6) +P2^(a)^(6) + P3^(a) ^(6) ■ 



(47) 



Consequently, we can interpret e as the energy density of the halo and p\ , P2 and P3 as the pressure in the principal 
directions of the halo. So, we have that the energy density of the halo is given by 



e = M 



(0)(0) 



whereas for the principal pressure we have 



P = Pi = P2 = ~P3 = M (l)(l), 



and 



<P> = 



P1+P2+ P3 



P 

3 



(48) 



(49) 



(50) 



is the average value of the pressure. 

Accordingly, we have an anisotropic fluid with a non-barotropic equation of state, which can be written as 



P = P(f, e) 



f-h 



3/ 



(51) 



in such a way that the pressure not only depends on the energy density but also on the gravitational and electromag- 
netic fields through the function /. Now, it is easy to see that, in order that the material energy-momentum tensor 
satisfy all the energy conditions, the function U must be constrained by 



U< 



(2 - fc 7 )fc 3 fc 4 
k 7 



(52) 



where we used the constrain on k$ in Eq. (|38l) needed to have an asymptotically flat spacctimc at infinity. 
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B. The Kuzmin-like solution 



One can find many different models for a relativistic thin disk surrounded by a material electromagnetizcd halo 
by choosing different kind of solutions U of the Laplace equation. Let us consider the particular case of Kuzmin's 
solution [Ej, [55j 

777 

U = - r= W , (a,m>0). (53) 

y/r 1 + (\z\ + a) z 

At points with z < 0, U is identical to the potential of a point mass m located at the point (r, z) = (0, —a), and when 
z > 0, f/ coincides with the potential generated by a point mass at (0, a). Hence V 2 C/ must vanish everywhere except 
on the plane z = 0. By applying Gauss theorem to a flat volume that contains a small portion of the plane z = 0, we 
conclude that U is generated by the surface density of a Newtonian mass 

am 

pM = 0)= 27r(r2 + Q2)3/2 - (54) 

Furthermore, without loss of generality we can choose the constants fcs = ki = 1. Then, for the metric potential we 
have 



e « _ ^ + <|i| + i)- , (55) 

y/f 2 + p| + l) 2 - m 

where we introduced the dimensionlcss variables r = r/a, z = z/a, and m = m/a. With this metric potential it 

is straightforward to calculate the dimensionless radial and axial components of the electric field, E r = aE r and 
E z = a 2 E z , and the dimensionless radial and axial components of the magnetic field, B r = B r and B z = B z . We 
obtain 

E r = f 7, (56a) 

k{f 2 + {\~z\ + 1 ) 2 } {v/f 2 + + l) 2 - mj 

£, = - — klM f [ i t l \. — r , (56b) 

/cz - 



i {f 2 + (|z| + f ) 2 } { v/r 2 + (|z| + l) 2 - m} ' 



B r = (56c) 

£;{f 2 + (|z| + l) 2 } 3/2 

4 = + 3/2 . (56d) 

fcz{f 2 + (|z| + l) 2 } 3/2 

were we have used the definitions E r = Ao, r , E z = Aq, z , B r = A^ z and B z = — A >r . 

Substituting fl53]) into gT]) and (l42bl) . we obtain for the dimensionless energy and charge density on surface of the 
disk, 6 == (26, o = aa, respectively, 

c(f) = ; 2m -, (57a) 

~/~\ 2klTfl fKVU\ 

<j(r) = ■ . (57b) 

v ' fc(f 2 + l)(\/FTl-m) v ' 

The substitution of (|55|) into P5]l and (Ti9")) allows us to write the dimensionless energy of the halo e(f, z) = a 2 e(r, z) 
and the dimensionless pressure of the halo p(r, z) = a 2 p{r, z) as 



mi) = t\ ,y^L-U, (58a) 



^ y/r^ + dzl + l) 2 1 

p(r ' zHr i ^ 2 + Pi + i) 2 -^ ' fc T (58b) 



where 



: {f 2 + (|z| + f ) 2 } {v/f 2 + (|5| + l) 2 - mV 
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(a) (b) 

FIG. 1. Dimensionless surface energy e and charge a densities as a function of r. In each case, we plot e(f) and a(f) for different 
values of the parameter in. First, we take rh = 0.45 (the bottom curve in each plot) and then 0.55, 0.75 and rh = 0.85 (the 
top curve in each plot). 



To fix the values of the constants that enter the solution, we consider the energy conditions at infinity, r , z — > oo, 
and at the center of the body, f, z — > 0. Then, it can be shown that all the energy conditions are satisfied at these 
extreme regions if 

m < 1 (59) 

for the disk and, additionally, 

- \Jk\ + l < k < yjkf + l (that is k 7 < 2) (60) 

for the halo. 

In Fig. []Ja), we show the dimensionless surface energy density on the disk e as a function of f and for different values 
of the parameter rh. First, we take rh = 0.45 (the bottom curve in the plot) and then 0.55, 0.75 and rh = 0.85 (the 
top curve in the plot). It can be seen that the energy density is everywhere positive fulfilling the energy conditions. 
It can be observed that for all the values of rh the maximum of the energy density occurs at the center of the disk and 
that it vanishes sufficiently fast as r increases. It can also be observed that the energy density in the central region 
of the disk increases as the values of the parameter rh increase. We have also plotted in Fig. [Tf b) the charge density 
a as a function of f. In each case, we plot a for different values of the parameter rh. We observe that the electric 
charge density has a behavior similar to that of the energy. This is consistent with the fact that the mass and the 
charge are more densely concentrated in the center of the disk. We also computed these functions for other values of 
the parameters rh in the interval (0, 1) and in all cases we found a similar behavior. 

In Fig. Ufa) and in Fig. [DJb), we illustrate the behavior of the surfaces and level curves of the matter density in 
the halo around of the disk (r > 0, z > 0) for the parameters rh = 075 and kj = 1. We can see that the energy is 
everywhere positive, its maximum occurs around the center of the disk, and it vanishes sufficiently fast as r increases. 

Notice that in the limit m — > 0, the gravitational and electromagnetic fields vanish identically, and the metric 
becomes flat. This is an important limiting case because it indicates that m determines the mass of the disk and the 
halo, and that the electromagnetic field exists only in connection with the disk-halo configuration. Moreover, when 
ki = 0, we obtain Aq — = const., a "purely magnetic" solution. Then, k± determines the electric charge of the 
distribution whereas kg = 1/k must be associated with the magnetic field. 

The energy condition rh < 1 or, equivalently, m < a imposes a maximum on the value of the mass parameter. 
Recall that in the Kuzmin solution a represents the distance along the axis between the equatorial plane and the 
point where the mass m is situated; accordingly, a can be interpreted as a parameter determining a proper length 
for the configuration. Then, the inequality m/a < 1 represents a condition on the "specific mass" of the system. 
This resembles the well-known Chandrasekhar limit for a spherically symmetric mass distribution, stating that the 
condition mass/radius < 4/9 must be satisfied in order to avoid gravitational collapse. The energy condition for the 
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FIG. 2. Surface plot an level curves of the energy density i on the exterior halo as a function of f and z with parameters 
rh = 0.75 and kj = 1. 




FIG. 3. Surface plot and level curves of the radial pressure p on the exterior halo as a function of f and z with parameters 
rh — 0.75 and k-j — 1. 

halo k 2 < 1 + k\ also represents a relationship between the parameters that characterize the electric and magnetic 
properties of the system. 

We conclude that all the relevant quantities show a physically reasonable behavior within the allowed range of 
values of the parameters. This indicates that the solution presented here can be used to describe the gravitational 
field of a static thin disk surrounded by a material halo with a non-trivial electromagnetic field. 



C. Singular behavior of the Kuzmin-like solution 

In order to study the singularities that could be present in the space described by the Kuzmin-like solution derived 
in the last subsection, we compute the most important quadratic curvature scalars, namely, the Krctschmann /Cj, the 
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0.2 0.4 0.6 0.8 1 



(a) 



(b) 



FIG. 4. Surface plot an level curves of the average value of the radial pressure < p > on the exterior halo as a function of r 
and z with parameters rh — 0.75 and &7 = 1. 



Chern-Pontryagin Kn and the Euler invariants Km defined as [56j 

y jjabed p 

/Vr — -K ftabed, 

ab nijed J3 



v- T* placed p 

A- J/ — [ ilj -fta&cd — 



c abij ,cdkl p p 
IT r* p*laf>cd p fc fc fi-ijklfiabcd 
I^III = [ it I ilabcd — • 

As for the Maxwell field, we consider the electromagnetic invariants 

Tt = F ab F ab 

77 rp rp*ab 



(61a) 
(61b) 

(61c) 



(62a) 
(62b) 



Here g = det (<?„&), e abcd is the Levi-Civita symbol and the asterisk denotes the dual operation. By using the solution 
([55]) we can cast these invariants as 



Ki(r,z) = 



|48[(|5| + l) 2 +r 2 ] -32^^(151 + 1)2+^ + 11m 2 } 



4[(|5| + l) 2 +f 2 ] y/(\z\ + l) 2 + f 2 - m 



Kn(f,z) = 0, 



K in (f,z) 



16m 2 {3[(|2| + l) 2 + r 2 ] - 2rhy/(\z\ + l) 2 + f 2 } 



whereas the electromagnetic invariants are 

•F, = 



z| + l) 2 +r 2 ] v /(|z| + l) 2 +r 2 -m 
2(1 - fc 2 )m 2 a 2 



fc 2 + l) 2 + f 2 ] [VP + l) 2 + f 2 - m 
-4fcim 2 a 2 



fc 2 [(ifi + i) 2 + f 2 ] ypi + iy+f 



(64a) 
(64b) 
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We see that there exists a singularity at the surface determined by the equation 



(|2| + l) 2 + f 2 =m 2 , 



(65) 



where all the non-trivial invariants diverge. The singularity equation allows real solutions only for m > 1. On the 
other hand, in the previous subsection we showed that the energy condition for the disk-halo configuration implies that 
m < 1. In fact, we can see from the expressions for the matter and charge density of the halo that they both diverge 
at the singular surface. Such divergencies are usually associated with the collapse of the gravitational configuration. 
This is an interesting result because it shows that the spacetime becomes singular as soon as the energy condition is 
violated. From a physical point of view this means that a disk-halo system can be constructed only if the condition 
m/a < 1 is satisfied. As noticed in the last subsection, this condition implies an upper bound on the value of the 
mass parameter m which is determined by the position of the mass along the symmetry axis. This, again, might be 
interpreted as a Chandrasekhar-like limit for the disk-halo system. 



In this work, we derived a rclativistic model describing a thin disk surrounded by a halo in presence of an electro- 
magnetic field. The model was obtained by solving the Einstein-Maxwell equations on a particular conformastatic 
spacetime in which only one independent metric function appears. For the energy-momentum tensor we used the 
distributional approach, and represented it as the sum of two distributional contributions, one due to the electro- 
magnetic part and the other associated with a matter distribution. These assumptions allowed us to derive explicit 
expressions for the energy, pressure, electric current and electromagnetic field of the disk region and the halo as well. 
The main point of this approach is that it allows to write the gravitational and electromagnetic potentials in terms 
of a solution of Laplace's equation. 

As a particular example, we used one of the simplest solutions of Laplace's equation, known as the Kuzmin so- 
lution, which contains two independent parameters, namely, the mass m and the parameter a that determines the 
proper length of the mass system. The resulting Kuzmin-like solution contains four independent parameters which 
determine the mass, proper length of the seed Kuzmin solution, electric charge and the magnetic field. The solution 
is asymptotically flat in general and turns out to be free of singularities if the ratio m/a < 1, an inequality that also 
guarantees the fulfillment of all the energy conditions. We interpret this condition as a Chandrasekhar-like limit for 
the disk-halo system. 

Since all the relevant quantities show a physically reasonable behavior within the range m/a < 1, we conclude 
that the solution presented here can be used to describe the gravitational and electromagnetic fields of a think disk 
surrounded by a halo in the presence of an electromagnetic field. 

Furthermore, with the particular solution of Laplace's equation here choosen, the matter distribution both of the 
disk and the halo vanishes sufficiently fast as r and z increase. This can be considered as indicating that the sources, 
the halo and the disk, are of finite size. Now, this reasonable behavior is a consequence of the behavior of the Kuzmin 
solution, which describes a gravitational potential that have a negative maximum at the center of the source and 
then decreases in absolute value, vanishing at infinity. Accordingly, we can expect that in a more general case a 
similar behavior will be obtained if solutions of the Laplace equation are considered that correspond to concentrated 
Newtonian sources with decreasing gravitational potential and gravitational field. 
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Appendix A: Appendix 



For the general comformastatic metric 



ds 2 = -c^dt 2 + e 2 ^[r 2 dip 2 + dr 2 + dz 2 



(Al) 
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the nonzero components of the Einstein tensor read 



Gto = ~e 2 ^~^{2(^ rr + ^ zz + -W) + i>\ + ^ 2 J, (A2a) 

r 

Gil = r 2 {4>,rr + (f 2 r + <t>,zz + <f>% + ^,rr + lp,zz}, (A2b) 

G^2 = 4>,zz + -<f>, r + (f> 2 z + 1p,zz + -4>,r + + 2</>,rV>,r, (A2c) 

G 23 = "f" - <$>,A,z + <l>,r^,z + §,A,r ~ '4 \rz + i>,r^,z, (A2d) 

Gf 3 = 0, rr + -<t> >r + (f) 2 r + VVrr + -1p,r + 1p% + ?>4> A * ■ ( A2 °) 

Furthermore, the nonzero components of the electromagnetic tensor (|10[) are given by 

^0 = \c- 2 HAl, r + Al z + r-V^)(A 2 r + A% )}, (A3a) 

= e-^(A , r A, r + A , Z A, Z ), (A3b) 

E± = \e-^{A 2 r + A% + r 2 e- 2 ^ « r + A^)}, (A3c) 

Et 2 = \{-e- 24, {Al r A 2 Z ) + r- 2 e-^{A 2 r - A 2 Z )}, (A3d) 

E± = -e-^Aa^Ao^ + r' 2 'e -2 * 'A, r A, g , (A3e) 

^3 = \{e~ H {Al r - A 2 Q J r- 2 e- 2 \A 2 r - A 2 Z )} , (A3f) 

where, all the quantities are evaluated in the z > and z < regions. Consequently, from the general equation 

G± b = Rtt ^.9a 6 i? ± = E% + M±, (A4) 
and Eqs. (|A2[) and (|A3[) . we have for the nonzero components of the matter energy- momentum tensor 

= -e 2( ^-^ ) {2VV + VV- VV> + ^e-^VAo- VA a + ^r~ 2 e _a *V A- VA}, (A5a) 

= -e-^VAo- VA, (A5b) 

= r 2 {V 2 {4> + V) - ~{4> + i/>),r + V0- - ^e^VAo- WA Q - -r~ 2 e'^W A- VA}, (A5c) 

M± = V 2 (0 + V) " (0 + V),rr + 2 z + ^ + 20,r^r + \e~ 2 *{A 2 . r - A\ z ) - ^ er 2 ^\A 2 r - A 2 Z ), (A5d) 

Mm = ~0 + V0>« -4>A<t> ~ ^),« + VvW> + 1>l* + e- 2 ^A . r A^ z - r^e^A^A,, (A5e) 
M 3 1 = V 2 (^ + V) - (4> + V0,« + 2 ,. + V,* + W,ztz - \e- 2 *{Al r - A 2 J + l -r- 2 e- 2 ^(A 2 r - A 2 Z ). (A5f) 



Alternatively, the components of the MEMT in the local orthonormal tetrad for the general metric (|A1[) can be 
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expressed as 



±r-*e-**{A* r -A%)}, 



M (o)(o) = -e-^W^rr + t** + \i\r) + ^ r + V> 2 + {A 2 ^ + A 2 , J + l - r - 2 e-^{A 2 r + A 2 Z , )}, (A6a) 

M± )(1) = - r -i e -(^) (A 0>r A tT + A 0>z A tZ ), (A6b) 

M± )(1) = e- 2 ^{0, rr + </> 2 r + + 4>l + trr + t*> ~ \f~ 2 ^Kr + ^o,J " \r~ 2 e^' {A 2 r + A%)}, (A6c) 

M (2)(2) = ^'R- + + <t>% + tzz + + V.r + 2^, r V,r + ^~ 2 Vo,r ~ ( A6d ) 
1 

— 1 

2 

M (2)(3) = e ~ 2,/ '{-^r Z - ^r^,* + 0, rV>,* + 0,,^,r " + ^, r ^,* + e^yV Aq,* - r^e -2 ^, r A J, (A6e) 

M (3)(3) = e_2V '{^,- + Vr + 4 + trr + U\r + tf z + 20,,V,« ~ ^"^(A 2 , r ~ A 2 J (A6f) 

+ \r- 2 e- 2 ^A 2 ,-A 2 z )} 1 

Jf 0) = re^{V 2 A -VAo-V(</> -</>)}, (A6g) 

0/1 

J ( ± = -e^{V 2 A + VA • V(</> - V) - — }• (A6h) 
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